3.6.67 \(\int \frac {(f+g x)^2}{(d+e x)^3 (d^2-e^2 x^2)^2} \, dx\) [567]

3.6.67.1 Optimal result
3.6.67.2 Mathematica [A] (verified)
3.6.67.3 Rubi [A] (verified)
3.6.67.4 Maple [A] (verified)
3.6.67.5 Fricas [B] (verification not implemented)
3.6.67.6 Sympy [B] (verification not implemented)
3.6.67.7 Maxima [A] (verification not implemented)
3.6.67.8 Giac [A] (verification not implemented)
3.6.67.9 Mupad [B] (verification not implemented)

3.6.67.1 Optimal result

Integrand size = 29, antiderivative size = 178 \[ \int \frac {(f+g x)^2}{(d+e x)^3 \left (d^2-e^2 x^2\right )^2} \, dx=\frac {(e f+d g)^2}{32 d^5 e^3 (d-e x)}-\frac {(e f-d g)^2}{16 d^2 e^3 (d+e x)^4}-\frac {e^2 f^2-d^2 g^2}{12 d^3 e^3 (d+e x)^3}-\frac {(3 e f-d g) (e f+d g)}{32 d^4 e^3 (d+e x)^2}-\frac {f (e f+d g)}{8 d^5 e^2 (d+e x)}+\frac {(e f+d g) (5 e f+d g) \text {arctanh}\left (\frac {e x}{d}\right )}{32 d^6 e^3} \]

output
1/32*(d*g+e*f)^2/d^5/e^3/(-e*x+d)-1/16*(-d*g+e*f)^2/d^2/e^3/(e*x+d)^4+1/12 
*(d^2*g^2-e^2*f^2)/d^3/e^3/(e*x+d)^3-1/32*(-d*g+3*e*f)*(d*g+e*f)/d^4/e^3/( 
e*x+d)^2-1/8*f*(d*g+e*f)/d^5/e^2/(e*x+d)+1/32*(d*g+e*f)*(d*g+5*e*f)*arctan 
h(e*x/d)/d^6/e^3
 
3.6.67.2 Mathematica [A] (verified)

Time = 0.09 (sec) , antiderivative size = 195, normalized size of antiderivative = 1.10 \[ \int \frac {(f+g x)^2}{(d+e x)^3 \left (d^2-e^2 x^2\right )^2} \, dx=\frac {\frac {6 d (e f+d g)^2}{d-e x}-\frac {12 d^4 (e f-d g)^2}{(d+e x)^4}+\frac {16 d^3 \left (-e^2 f^2+d^2 g^2\right )}{(d+e x)^3}+\frac {6 d^2 \left (-3 e^2 f^2-2 d e f g+d^2 g^2\right )}{(d+e x)^2}-\frac {24 d e f (e f+d g)}{d+e x}-3 \left (5 e^2 f^2+6 d e f g+d^2 g^2\right ) \log (d-e x)+3 \left (5 e^2 f^2+6 d e f g+d^2 g^2\right ) \log (d+e x)}{192 d^6 e^3} \]

input
Integrate[(f + g*x)^2/((d + e*x)^3*(d^2 - e^2*x^2)^2),x]
 
output
((6*d*(e*f + d*g)^2)/(d - e*x) - (12*d^4*(e*f - d*g)^2)/(d + e*x)^4 + (16* 
d^3*(-(e^2*f^2) + d^2*g^2))/(d + e*x)^3 + (6*d^2*(-3*e^2*f^2 - 2*d*e*f*g + 
 d^2*g^2))/(d + e*x)^2 - (24*d*e*f*(e*f + d*g))/(d + e*x) - 3*(5*e^2*f^2 + 
 6*d*e*f*g + d^2*g^2)*Log[d - e*x] + 3*(5*e^2*f^2 + 6*d*e*f*g + d^2*g^2)*L 
og[d + e*x])/(192*d^6*e^3)
 
3.6.67.3 Rubi [A] (verified)

Time = 0.40 (sec) , antiderivative size = 178, normalized size of antiderivative = 1.00, number of steps used = 3, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.103, Rules used = {639, 99, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {(f+g x)^2}{(d+e x)^3 \left (d^2-e^2 x^2\right )^2} \, dx\)

\(\Big \downarrow \) 639

\(\displaystyle \int \frac {(f+g x)^2}{(d-e x)^2 (d+e x)^5}dx\)

\(\Big \downarrow \) 99

\(\displaystyle \int \left (\frac {(d g+e f)^2}{32 d^5 e^2 (d-e x)^2}+\frac {f (d g+e f)}{8 d^5 e (d+e x)^2}+\frac {(3 e f-d g) (d g+e f)}{16 d^4 e^2 (d+e x)^3}+\frac {(d g-e f)^2}{4 d^2 e^2 (d+e x)^5}+\frac {(d g+e f) (d g+5 e f)}{32 d^5 e^2 \left (d^2-e^2 x^2\right )}+\frac {e^2 f^2-d^2 g^2}{4 d^3 e^2 (d+e x)^4}\right )dx\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {\text {arctanh}\left (\frac {e x}{d}\right ) (d g+e f) (d g+5 e f)}{32 d^6 e^3}+\frac {(d g+e f)^2}{32 d^5 e^3 (d-e x)}-\frac {f (d g+e f)}{8 d^5 e^2 (d+e x)}-\frac {(3 e f-d g) (d g+e f)}{32 d^4 e^3 (d+e x)^2}-\frac {(e f-d g)^2}{16 d^2 e^3 (d+e x)^4}-\frac {e^2 f^2-d^2 g^2}{12 d^3 e^3 (d+e x)^3}\)

input
Int[(f + g*x)^2/((d + e*x)^3*(d^2 - e^2*x^2)^2),x]
 
output
(e*f + d*g)^2/(32*d^5*e^3*(d - e*x)) - (e*f - d*g)^2/(16*d^2*e^3*(d + e*x) 
^4) - (e^2*f^2 - d^2*g^2)/(12*d^3*e^3*(d + e*x)^3) - ((3*e*f - d*g)*(e*f + 
 d*g))/(32*d^4*e^3*(d + e*x)^2) - (f*(e*f + d*g))/(8*d^5*e^2*(d + e*x)) + 
((e*f + d*g)*(5*e*f + d*g)*ArcTanh[(e*x)/d])/(32*d^6*e^3)
 

3.6.67.3.1 Defintions of rubi rules used

rule 99
Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_) 
)^(p_), x_] :> Int[ExpandIntegrand[(a + b*x)^m*(c + d*x)^n*(e + f*x)^p, x], 
 x] /; FreeQ[{a, b, c, d, e, f, p}, x] && IntegersQ[m, n] && (IntegerQ[p] | 
| (GtQ[m, 0] && GeQ[n, -1]))
 

rule 639
Int[((c_) + (d_.)*(x_))^(m_.)*((e_) + (f_.)*(x_))^(n_.)*((a_) + (b_.)*(x_)^ 
2)^(p_.), x_Symbol] :> Int[(c + d*x)^(m + p)*(e + f*x)^n*(a/c + (b/d)*x)^p, 
 x] /; FreeQ[{a, b, c, d, e, f, m, n, p}, x] && EqQ[b*c^2 + a*d^2, 0] && (I 
ntegerQ[p] || (GtQ[a, 0] && GtQ[c, 0] &&  !IntegerQ[m]))
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 
3.6.67.4 Maple [A] (verified)

Time = 0.42 (sec) , antiderivative size = 241, normalized size of antiderivative = 1.35

method result size
default \(\frac {\left (-d^{2} g^{2}-6 d e f g -5 e^{2} f^{2}\right ) \ln \left (-e x +d \right )}{64 e^{3} d^{6}}+\frac {d^{2} g^{2}+2 d e f g +e^{2} f^{2}}{32 e^{3} d^{5} \left (-e x +d \right )}+\frac {\left (d^{2} g^{2}+6 d e f g +5 e^{2} f^{2}\right ) \ln \left (e x +d \right )}{64 e^{3} d^{6}}-\frac {-d^{2} g^{2}+e^{2} f^{2}}{12 e^{3} d^{3} \left (e x +d \right )^{3}}-\frac {-d^{2} g^{2}+2 d e f g +3 e^{2} f^{2}}{32 e^{3} d^{4} \left (e x +d \right )^{2}}-\frac {d^{2} g^{2}-2 d e f g +e^{2} f^{2}}{16 d^{2} e^{3} \left (e x +d \right )^{4}}-\frac {f \left (d g +e f \right )}{8 d^{5} e^{2} \left (e x +d \right )}\) \(241\)
norman \(\frac {\frac {\left (25 d^{2} g^{2}+54 d e f g -19 e^{2} f^{2}\right ) x^{3}}{96 d^{4}}-\frac {\left (3 d^{2} g^{2}-14 d e f g -33 e^{2} f^{2}\right ) x^{2}}{32 e \,d^{3}}+\frac {3 e \left (3 d^{2} g^{2}+2 d e f g -9 e^{2} f^{2}\right ) x^{4}}{32 d^{5}}+\frac {e^{2} \left (d^{2} g^{2}-4 e^{2} f^{2}\right ) x^{5}}{12 d^{6}}-\frac {\left (d^{2} g^{2}+6 d e f g -27 e^{2} f^{2}\right ) x}{32 d^{2} e^{2}}}{\left (e x +d \right )^{4} \left (-e x +d \right )}-\frac {\left (d^{2} g^{2}+6 d e f g +5 e^{2} f^{2}\right ) \ln \left (-e x +d \right )}{64 e^{3} d^{6}}+\frac {\left (d^{2} g^{2}+6 d e f g +5 e^{2} f^{2}\right ) \ln \left (e x +d \right )}{64 e^{3} d^{6}}\) \(247\)
risch \(\frac {\frac {e \left (d^{2} g^{2}+6 d e f g +5 e^{2} f^{2}\right ) x^{4}}{32 d^{5}}+\frac {3 \left (d^{2} g^{2}+6 d e f g +5 e^{2} f^{2}\right ) x^{3}}{32 d^{4}}+\frac {7 \left (d^{2} g^{2}+6 d e f g +5 e^{2} f^{2}\right ) x^{2}}{96 d^{3} e}+\frac {\left (7 d^{2} g^{2}-6 d e f g -5 e^{2} f^{2}\right ) x}{32 d^{2} e^{2}}+\frac {d^{2} g^{2}-4 e^{2} f^{2}}{12 d \,e^{3}}}{\left (e x +d \right )^{3} \left (-e^{2} x^{2}+d^{2}\right )}-\frac {\ln \left (-e x +d \right ) g^{2}}{64 e^{3} d^{4}}-\frac {3 \ln \left (-e x +d \right ) f g}{32 e^{2} d^{5}}-\frac {5 \ln \left (-e x +d \right ) f^{2}}{64 e \,d^{6}}+\frac {\ln \left (e x +d \right ) g^{2}}{64 e^{3} d^{4}}+\frac {3 \ln \left (e x +d \right ) f g}{32 e^{2} d^{5}}+\frac {5 \ln \left (e x +d \right ) f^{2}}{64 e \,d^{6}}\) \(278\)
parallelrisch \(-\frac {-15 \ln \left (e x +d \right ) x^{5} e^{7} f^{2}+15 \ln \left (e x -d \right ) x^{5} e^{7} f^{2}+162 x \,d^{4} e^{3} f^{2}+54 \ln \left (e x +d \right ) x \,d^{5} e^{2} f g +3 \ln \left (e x -d \right ) x^{5} d^{2} e^{5} g^{2}-3 \ln \left (e x +d \right ) x^{5} d^{2} e^{5} g^{2}+9 \ln \left (e x -d \right ) x^{4} d^{3} e^{4} g^{2}+45 \ln \left (e x -d \right ) x^{4} d \,e^{6} f^{2}-9 \ln \left (e x +d \right ) x^{4} d^{3} e^{4} g^{2}-45 \ln \left (e x +d \right ) x^{4} d \,e^{6} f^{2}+6 \ln \left (e x -d \right ) x^{3} d^{4} e^{3} g^{2}+30 \ln \left (e x -d \right ) x^{3} d^{2} e^{5} f^{2}-9 \ln \left (e x -d \right ) x \,d^{6} e \,g^{2}-45 \ln \left (e x -d \right ) x \,d^{4} e^{3} f^{2}-18 \ln \left (e x -d \right ) d^{6} e f g -6 x \,d^{6} e \,g^{2}+198 x^{2} d^{3} e^{4} f^{2}+50 x^{3} d^{4} e^{3} g^{2}-38 x^{3} d^{2} e^{5} f^{2}+54 x^{4} d^{3} e^{4} g^{2}-162 x^{4} d \,e^{6} f^{2}+16 x^{5} d^{2} e^{5} g^{2}-15 \ln \left (e x -d \right ) d^{5} e^{2} f^{2}-18 x^{2} d^{5} e^{2} g^{2}-54 \ln \left (e x -d \right ) x \,d^{5} e^{2} f g +15 \ln \left (e x +d \right ) d^{5} e^{2} f^{2}-6 \ln \left (e x +d \right ) x^{3} d^{4} e^{3} g^{2}-30 \ln \left (e x +d \right ) x^{3} d^{2} e^{5} f^{2}-6 \ln \left (e x -d \right ) x^{2} d^{5} e^{2} g^{2}-30 \ln \left (e x -d \right ) x^{2} d^{3} e^{4} f^{2}+6 \ln \left (e x +d \right ) x^{2} d^{5} e^{2} g^{2}+30 \ln \left (e x +d \right ) x^{2} d^{3} e^{4} f^{2}+9 \ln \left (e x +d \right ) x \,d^{6} e \,g^{2}+45 \ln \left (e x +d \right ) x \,d^{4} e^{3} f^{2}+18 \ln \left (e x +d \right ) d^{6} e f g -36 x \,d^{5} e^{2} f g -3 \ln \left (e x -d \right ) d^{7} g^{2}+18 \ln \left (e x -d \right ) x^{5} d \,e^{6} f g -18 \ln \left (e x +d \right ) x^{5} d \,e^{6} f g -64 x^{5} e^{7} f^{2}+3 \ln \left (e x +d \right ) d^{7} g^{2}+84 x^{2} d^{4} e^{3} f g +108 x^{3} d^{3} e^{4} f g +36 x^{4} d^{2} e^{5} f g +54 \ln \left (e x -d \right ) x^{4} d^{2} e^{5} f g -54 \ln \left (e x +d \right ) x^{4} d^{2} e^{5} f g +36 \ln \left (e x -d \right ) x^{3} d^{3} e^{4} f g -36 \ln \left (e x +d \right ) x^{3} d^{3} e^{4} f g -36 \ln \left (e x -d \right ) x^{2} d^{4} e^{3} f g +36 \ln \left (e x +d \right ) x^{2} d^{4} e^{3} f g}{192 e^{3} d^{6} \left (e^{2} x^{2}-d^{2}\right ) \left (e x +d \right )^{3}}\) \(897\)

input
int((g*x+f)^2/(e*x+d)^3/(-e^2*x^2+d^2)^2,x,method=_RETURNVERBOSE)
 
output
1/64*(-d^2*g^2-6*d*e*f*g-5*e^2*f^2)/e^3/d^6*ln(-e*x+d)+1/32*(d^2*g^2+2*d*e 
*f*g+e^2*f^2)/e^3/d^5/(-e*x+d)+1/64/e^3*(d^2*g^2+6*d*e*f*g+5*e^2*f^2)/d^6* 
ln(e*x+d)-1/12*(-d^2*g^2+e^2*f^2)/e^3/d^3/(e*x+d)^3-1/32*(-d^2*g^2+2*d*e*f 
*g+3*e^2*f^2)/e^3/d^4/(e*x+d)^2-1/16*(d^2*g^2-2*d*e*f*g+e^2*f^2)/d^2/e^3/( 
e*x+d)^4-1/8*f*(d*g+e*f)/d^5/e^2/(e*x+d)
 
3.6.67.5 Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 648 vs. \(2 (167) = 334\).

Time = 0.35 (sec) , antiderivative size = 648, normalized size of antiderivative = 3.64 \[ \int \frac {(f+g x)^2}{(d+e x)^3 \left (d^2-e^2 x^2\right )^2} \, dx=\frac {64 \, d^{5} e^{2} f^{2} - 16 \, d^{7} g^{2} - 6 \, {\left (5 \, d e^{6} f^{2} + 6 \, d^{2} e^{5} f g + d^{3} e^{4} g^{2}\right )} x^{4} - 18 \, {\left (5 \, d^{2} e^{5} f^{2} + 6 \, d^{3} e^{4} f g + d^{4} e^{3} g^{2}\right )} x^{3} - 14 \, {\left (5 \, d^{3} e^{4} f^{2} + 6 \, d^{4} e^{3} f g + d^{5} e^{2} g^{2}\right )} x^{2} + 6 \, {\left (5 \, d^{4} e^{3} f^{2} + 6 \, d^{5} e^{2} f g - 7 \, d^{6} e g^{2}\right )} x - 3 \, {\left (5 \, d^{5} e^{2} f^{2} + 6 \, d^{6} e f g + d^{7} g^{2} - {\left (5 \, e^{7} f^{2} + 6 \, d e^{6} f g + d^{2} e^{5} g^{2}\right )} x^{5} - 3 \, {\left (5 \, d e^{6} f^{2} + 6 \, d^{2} e^{5} f g + d^{3} e^{4} g^{2}\right )} x^{4} - 2 \, {\left (5 \, d^{2} e^{5} f^{2} + 6 \, d^{3} e^{4} f g + d^{4} e^{3} g^{2}\right )} x^{3} + 2 \, {\left (5 \, d^{3} e^{4} f^{2} + 6 \, d^{4} e^{3} f g + d^{5} e^{2} g^{2}\right )} x^{2} + 3 \, {\left (5 \, d^{4} e^{3} f^{2} + 6 \, d^{5} e^{2} f g + d^{6} e g^{2}\right )} x\right )} \log \left (e x + d\right ) + 3 \, {\left (5 \, d^{5} e^{2} f^{2} + 6 \, d^{6} e f g + d^{7} g^{2} - {\left (5 \, e^{7} f^{2} + 6 \, d e^{6} f g + d^{2} e^{5} g^{2}\right )} x^{5} - 3 \, {\left (5 \, d e^{6} f^{2} + 6 \, d^{2} e^{5} f g + d^{3} e^{4} g^{2}\right )} x^{4} - 2 \, {\left (5 \, d^{2} e^{5} f^{2} + 6 \, d^{3} e^{4} f g + d^{4} e^{3} g^{2}\right )} x^{3} + 2 \, {\left (5 \, d^{3} e^{4} f^{2} + 6 \, d^{4} e^{3} f g + d^{5} e^{2} g^{2}\right )} x^{2} + 3 \, {\left (5 \, d^{4} e^{3} f^{2} + 6 \, d^{5} e^{2} f g + d^{6} e g^{2}\right )} x\right )} \log \left (e x - d\right )}{192 \, {\left (d^{6} e^{8} x^{5} + 3 \, d^{7} e^{7} x^{4} + 2 \, d^{8} e^{6} x^{3} - 2 \, d^{9} e^{5} x^{2} - 3 \, d^{10} e^{4} x - d^{11} e^{3}\right )}} \]

input
integrate((g*x+f)^2/(e*x+d)^3/(-e^2*x^2+d^2)^2,x, algorithm="fricas")
 
output
1/192*(64*d^5*e^2*f^2 - 16*d^7*g^2 - 6*(5*d*e^6*f^2 + 6*d^2*e^5*f*g + d^3* 
e^4*g^2)*x^4 - 18*(5*d^2*e^5*f^2 + 6*d^3*e^4*f*g + d^4*e^3*g^2)*x^3 - 14*( 
5*d^3*e^4*f^2 + 6*d^4*e^3*f*g + d^5*e^2*g^2)*x^2 + 6*(5*d^4*e^3*f^2 + 6*d^ 
5*e^2*f*g - 7*d^6*e*g^2)*x - 3*(5*d^5*e^2*f^2 + 6*d^6*e*f*g + d^7*g^2 - (5 
*e^7*f^2 + 6*d*e^6*f*g + d^2*e^5*g^2)*x^5 - 3*(5*d*e^6*f^2 + 6*d^2*e^5*f*g 
 + d^3*e^4*g^2)*x^4 - 2*(5*d^2*e^5*f^2 + 6*d^3*e^4*f*g + d^4*e^3*g^2)*x^3 
+ 2*(5*d^3*e^4*f^2 + 6*d^4*e^3*f*g + d^5*e^2*g^2)*x^2 + 3*(5*d^4*e^3*f^2 + 
 6*d^5*e^2*f*g + d^6*e*g^2)*x)*log(e*x + d) + 3*(5*d^5*e^2*f^2 + 6*d^6*e*f 
*g + d^7*g^2 - (5*e^7*f^2 + 6*d*e^6*f*g + d^2*e^5*g^2)*x^5 - 3*(5*d*e^6*f^ 
2 + 6*d^2*e^5*f*g + d^3*e^4*g^2)*x^4 - 2*(5*d^2*e^5*f^2 + 6*d^3*e^4*f*g + 
d^4*e^3*g^2)*x^3 + 2*(5*d^3*e^4*f^2 + 6*d^4*e^3*f*g + d^5*e^2*g^2)*x^2 + 3 
*(5*d^4*e^3*f^2 + 6*d^5*e^2*f*g + d^6*e*g^2)*x)*log(e*x - d))/(d^6*e^8*x^5 
 + 3*d^7*e^7*x^4 + 2*d^8*e^6*x^3 - 2*d^9*e^5*x^2 - 3*d^10*e^4*x - d^11*e^3 
)
 
3.6.67.6 Sympy [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 376 vs. \(2 (162) = 324\).

Time = 0.80 (sec) , antiderivative size = 376, normalized size of antiderivative = 2.11 \[ \int \frac {(f+g x)^2}{(d+e x)^3 \left (d^2-e^2 x^2\right )^2} \, dx=\frac {- 8 d^{6} g^{2} + 32 d^{4} e^{2} f^{2} + x^{4} \left (- 3 d^{2} e^{4} g^{2} - 18 d e^{5} f g - 15 e^{6} f^{2}\right ) + x^{3} \left (- 9 d^{3} e^{3} g^{2} - 54 d^{2} e^{4} f g - 45 d e^{5} f^{2}\right ) + x^{2} \left (- 7 d^{4} e^{2} g^{2} - 42 d^{3} e^{3} f g - 35 d^{2} e^{4} f^{2}\right ) + x \left (- 21 d^{5} e g^{2} + 18 d^{4} e^{2} f g + 15 d^{3} e^{3} f^{2}\right )}{- 96 d^{10} e^{3} - 288 d^{9} e^{4} x - 192 d^{8} e^{5} x^{2} + 192 d^{7} e^{6} x^{3} + 288 d^{6} e^{7} x^{4} + 96 d^{5} e^{8} x^{5}} - \frac {\left (d g + e f\right ) \left (d g + 5 e f\right ) \log {\left (- \frac {d \left (d g + e f\right ) \left (d g + 5 e f\right )}{e \left (d^{2} g^{2} + 6 d e f g + 5 e^{2} f^{2}\right )} + x \right )}}{64 d^{6} e^{3}} + \frac {\left (d g + e f\right ) \left (d g + 5 e f\right ) \log {\left (\frac {d \left (d g + e f\right ) \left (d g + 5 e f\right )}{e \left (d^{2} g^{2} + 6 d e f g + 5 e^{2} f^{2}\right )} + x \right )}}{64 d^{6} e^{3}} \]

input
integrate((g*x+f)**2/(e*x+d)**3/(-e**2*x**2+d**2)**2,x)
 
output
(-8*d**6*g**2 + 32*d**4*e**2*f**2 + x**4*(-3*d**2*e**4*g**2 - 18*d*e**5*f* 
g - 15*e**6*f**2) + x**3*(-9*d**3*e**3*g**2 - 54*d**2*e**4*f*g - 45*d*e**5 
*f**2) + x**2*(-7*d**4*e**2*g**2 - 42*d**3*e**3*f*g - 35*d**2*e**4*f**2) + 
 x*(-21*d**5*e*g**2 + 18*d**4*e**2*f*g + 15*d**3*e**3*f**2))/(-96*d**10*e* 
*3 - 288*d**9*e**4*x - 192*d**8*e**5*x**2 + 192*d**7*e**6*x**3 + 288*d**6* 
e**7*x**4 + 96*d**5*e**8*x**5) - (d*g + e*f)*(d*g + 5*e*f)*log(-d*(d*g + e 
*f)*(d*g + 5*e*f)/(e*(d**2*g**2 + 6*d*e*f*g + 5*e**2*f**2)) + x)/(64*d**6* 
e**3) + (d*g + e*f)*(d*g + 5*e*f)*log(d*(d*g + e*f)*(d*g + 5*e*f)/(e*(d**2 
*g**2 + 6*d*e*f*g + 5*e**2*f**2)) + x)/(64*d**6*e**3)
 
3.6.67.7 Maxima [A] (verification not implemented)

Time = 0.20 (sec) , antiderivative size = 298, normalized size of antiderivative = 1.67 \[ \int \frac {(f+g x)^2}{(d+e x)^3 \left (d^2-e^2 x^2\right )^2} \, dx=\frac {32 \, d^{4} e^{2} f^{2} - 8 \, d^{6} g^{2} - 3 \, {\left (5 \, e^{6} f^{2} + 6 \, d e^{5} f g + d^{2} e^{4} g^{2}\right )} x^{4} - 9 \, {\left (5 \, d e^{5} f^{2} + 6 \, d^{2} e^{4} f g + d^{3} e^{3} g^{2}\right )} x^{3} - 7 \, {\left (5 \, d^{2} e^{4} f^{2} + 6 \, d^{3} e^{3} f g + d^{4} e^{2} g^{2}\right )} x^{2} + 3 \, {\left (5 \, d^{3} e^{3} f^{2} + 6 \, d^{4} e^{2} f g - 7 \, d^{5} e g^{2}\right )} x}{96 \, {\left (d^{5} e^{8} x^{5} + 3 \, d^{6} e^{7} x^{4} + 2 \, d^{7} e^{6} x^{3} - 2 \, d^{8} e^{5} x^{2} - 3 \, d^{9} e^{4} x - d^{10} e^{3}\right )}} + \frac {{\left (5 \, e^{2} f^{2} + 6 \, d e f g + d^{2} g^{2}\right )} \log \left (e x + d\right )}{64 \, d^{6} e^{3}} - \frac {{\left (5 \, e^{2} f^{2} + 6 \, d e f g + d^{2} g^{2}\right )} \log \left (e x - d\right )}{64 \, d^{6} e^{3}} \]

input
integrate((g*x+f)^2/(e*x+d)^3/(-e^2*x^2+d^2)^2,x, algorithm="maxima")
 
output
1/96*(32*d^4*e^2*f^2 - 8*d^6*g^2 - 3*(5*e^6*f^2 + 6*d*e^5*f*g + d^2*e^4*g^ 
2)*x^4 - 9*(5*d*e^5*f^2 + 6*d^2*e^4*f*g + d^3*e^3*g^2)*x^3 - 7*(5*d^2*e^4* 
f^2 + 6*d^3*e^3*f*g + d^4*e^2*g^2)*x^2 + 3*(5*d^3*e^3*f^2 + 6*d^4*e^2*f*g 
- 7*d^5*e*g^2)*x)/(d^5*e^8*x^5 + 3*d^6*e^7*x^4 + 2*d^7*e^6*x^3 - 2*d^8*e^5 
*x^2 - 3*d^9*e^4*x - d^10*e^3) + 1/64*(5*e^2*f^2 + 6*d*e*f*g + d^2*g^2)*lo 
g(e*x + d)/(d^6*e^3) - 1/64*(5*e^2*f^2 + 6*d*e*f*g + d^2*g^2)*log(e*x - d) 
/(d^6*e^3)
 
3.6.67.8 Giac [A] (verification not implemented)

Time = 0.28 (sec) , antiderivative size = 264, normalized size of antiderivative = 1.48 \[ \int \frac {(f+g x)^2}{(d+e x)^3 \left (d^2-e^2 x^2\right )^2} \, dx=\frac {{\left (5 \, e^{2} f^{2} + 6 \, d e f g + d^{2} g^{2}\right )} \log \left ({\left | e x + d \right |}\right )}{64 \, d^{6} e^{3}} - \frac {{\left (5 \, e^{2} f^{2} + 6 \, d e f g + d^{2} g^{2}\right )} \log \left ({\left | e x - d \right |}\right )}{64 \, d^{6} e^{3}} + \frac {32 \, d^{5} e^{2} f^{2} - 8 \, d^{7} g^{2} - 3 \, {\left (5 \, d e^{6} f^{2} + 6 \, d^{2} e^{5} f g + d^{3} e^{4} g^{2}\right )} x^{4} - 9 \, {\left (5 \, d^{2} e^{5} f^{2} + 6 \, d^{3} e^{4} f g + d^{4} e^{3} g^{2}\right )} x^{3} - 7 \, {\left (5 \, d^{3} e^{4} f^{2} + 6 \, d^{4} e^{3} f g + d^{5} e^{2} g^{2}\right )} x^{2} + 3 \, {\left (5 \, d^{4} e^{3} f^{2} + 6 \, d^{5} e^{2} f g - 7 \, d^{6} e g^{2}\right )} x}{96 \, {\left (e x + d\right )}^{4} {\left (e x - d\right )} d^{6} e^{3}} \]

input
integrate((g*x+f)^2/(e*x+d)^3/(-e^2*x^2+d^2)^2,x, algorithm="giac")
 
output
1/64*(5*e^2*f^2 + 6*d*e*f*g + d^2*g^2)*log(abs(e*x + d))/(d^6*e^3) - 1/64* 
(5*e^2*f^2 + 6*d*e*f*g + d^2*g^2)*log(abs(e*x - d))/(d^6*e^3) + 1/96*(32*d 
^5*e^2*f^2 - 8*d^7*g^2 - 3*(5*d*e^6*f^2 + 6*d^2*e^5*f*g + d^3*e^4*g^2)*x^4 
 - 9*(5*d^2*e^5*f^2 + 6*d^3*e^4*f*g + d^4*e^3*g^2)*x^3 - 7*(5*d^3*e^4*f^2 
+ 6*d^4*e^3*f*g + d^5*e^2*g^2)*x^2 + 3*(5*d^4*e^3*f^2 + 6*d^5*e^2*f*g - 7* 
d^6*e*g^2)*x)/((e*x + d)^4*(e*x - d)*d^6*e^3)
 
3.6.67.9 Mupad [B] (verification not implemented)

Time = 11.80 (sec) , antiderivative size = 274, normalized size of antiderivative = 1.54 \[ \int \frac {(f+g x)^2}{(d+e x)^3 \left (d^2-e^2 x^2\right )^2} \, dx=\frac {\frac {d^2\,g^2-4\,e^2\,f^2}{12\,d\,e^3}+\frac {3\,x^3\,\left (d^2\,g^2+6\,d\,e\,f\,g+5\,e^2\,f^2\right )}{32\,d^4}+\frac {e\,x^4\,\left (d^2\,g^2+6\,d\,e\,f\,g+5\,e^2\,f^2\right )}{32\,d^5}-\frac {x\,\left (-7\,d^2\,g^2+6\,d\,e\,f\,g+5\,e^2\,f^2\right )}{32\,d^2\,e^2}+\frac {7\,x^2\,\left (d^2\,g^2+6\,d\,e\,f\,g+5\,e^2\,f^2\right )}{96\,d^3\,e}}{d^5+3\,d^4\,e\,x+2\,d^3\,e^2\,x^2-2\,d^2\,e^3\,x^3-3\,d\,e^4\,x^4-e^5\,x^5}+\frac {\mathrm {atanh}\left (\frac {e\,x\,\left (d\,g+e\,f\right )\,\left (d\,g+5\,e\,f\right )}{d\,\left (d^2\,g^2+6\,d\,e\,f\,g+5\,e^2\,f^2\right )}\right )\,\left (d\,g+e\,f\right )\,\left (d\,g+5\,e\,f\right )}{32\,d^6\,e^3} \]

input
int((f + g*x)^2/((d^2 - e^2*x^2)^2*(d + e*x)^3),x)
 
output
((d^2*g^2 - 4*e^2*f^2)/(12*d*e^3) + (3*x^3*(d^2*g^2 + 5*e^2*f^2 + 6*d*e*f* 
g))/(32*d^4) + (e*x^4*(d^2*g^2 + 5*e^2*f^2 + 6*d*e*f*g))/(32*d^5) - (x*(5* 
e^2*f^2 - 7*d^2*g^2 + 6*d*e*f*g))/(32*d^2*e^2) + (7*x^2*(d^2*g^2 + 5*e^2*f 
^2 + 6*d*e*f*g))/(96*d^3*e))/(d^5 - e^5*x^5 - 3*d*e^4*x^4 + 2*d^3*e^2*x^2 
- 2*d^2*e^3*x^3 + 3*d^4*e*x) + (atanh((e*x*(d*g + e*f)*(d*g + 5*e*f))/(d*( 
d^2*g^2 + 5*e^2*f^2 + 6*d*e*f*g)))*(d*g + e*f)*(d*g + 5*e*f))/(32*d^6*e^3)